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ABSTRACT 


The phenomenon of mutual coupling in antenna arrays is a potential source of 
performance degradation, particularly in a highly congested environment such as the 
unmanned aerial vehicle where multitudes of arrays are accommodated in a constrained 
platform. This thesis investigates the effects of mutual coupling between elements in an 
antenna system via software simulations that use the Methods of Moments, a numerical 
technique that accounts for mutual coupling. The antenna configurations investigated 
comprise a single dipole and a three-dipole array, both above a simulated infinite perfect 
ground plane. 

The simulation results show that the active impedance of any element in an array 
depends on both the self and mutual impedances and that the terminal current phases, 
hence directivity, of an array are affected by mutual coupling such that the array pattern 
deviates from that defined by classical theoretical approach which assumes identical 
element patterns. Results have been obtained for a single dipole and a three-dipole array 
above a simulated ground plane. The examination of a small array provides a useful 
environment in which to develop, optimize and evaluate the radiating elements, thus 
providing better understanding of the effects of mutual coupling and facilitating the 
design of corresponding compensation techniques in practical antenna arrays. 
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EXECUTIVE SUMMARY 


In the classical approach of antenna theories, mutual coupling between elements 
is usually ignored. However, depending on the application, errors due to mutual coupling 
can be significant. The phenomenon of mutual coupling in antenna arrays is a potential 
source of performance degradation, particularly in a highly congested environment such 
as the unmanned aerial vehicle where multitudes of arrays are accommodated in a 
constrained platform. 

This thesis investigates the effects of mutual coupling between elements in an 
antenna system via software simulations that use the Methods of Moments (MM), a 
numerical technique that accounts for mutual coupling, as opposed to the classical 
approach, which does not. The antenna configurations investigated comprise a single 
dipole and a three-dipole array, both above a simulated infinite perfect ground plane. 

The frequency used in the simulations is 300 MHz. The basic dipole geometry 
consisted of a quadrilateral of length Z = 0.45/1 and width w= 0.005A . This length is 
chosen so that the current on the entire dipole is in phase, like in most practical 
applications. The length-to-width ratio is 90, which closely approximates an infinitely 
thin dipole. 

The antenna analysis software packages, Patch and NEC-Win Professional, are 
used to determine the active impedance of a centrally voltage-excited dipole in free space 
and its mutual impedance with a perfect ground. Although Patch and NEC-Win 
Professional use the MM for antenna analysis, differences inherent in the coding and 
assumptions are expected. Slight variations were therefore made to the basic dipole 
geometry in NEC Win Professional. The dipole was placed parallel as well as 
perpendicular to the ground plane and its distance above the ground plane was varied. In 
the presence of an infinite ground plane, the driving point impedance of the dipole is 
found to be a function of terminal currents and mutual impedance. The mutual impedance 
is due to coupling between the dipole and the ground plane, and this effect diminishes as 
the dipole is moved away from the ground plane. 
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The effects of mutual coupling on the phase, directivity and active element 
patterns of a three-dipole array was investigated for dipole separations of 0.5/1 and 
0.375/1, with the array oriented horizontally with respect to the ground and placed at a 
height 0.25/1 above it. In the transmit scenario, each dipole is centrally excited by a 
complex voltage of unit amplitude according to the dipole position in the array when the 
array is scanned. In the receive scenario, an incoming plane wave of unit amplitude is 
used. Mutual coupling changes the current magnitude, phase and distribution on the 
dipole, thereby ‘distorting’ the radiation vectors, giving rise to peaks and nulls in 
directions that are different from classical theoretical approximations. 

The active element pattern depends on the position of the fed element in the array. 
The importance of the results on active element pattern comes from the fact that 
measurements of the active element patterns can give more directly the desired scan 
characteristics of an array under test and are therefore useful in predicting the scan 
performance of large phased arrays. Thus, active element patterns can be used to deduce 
the radiation properties of larger arrays. 

The presence of mutual coupling variations in a finite antenna array causes 
amplitude and phase errors, the extent of which depends on the geometry and 
configuration of the array. Depending on the applications, such errors can be significant 
and result in severe system performance degradation, and must therefore be compensated 
for accordingly. In light of this, the research has also discussed further works with regard 
to compensating for the mutual coupling effects in dipole arrays. Even if these techniques 
do not result in the desired end-state of implementation in operational arrays, they can 
still be useful in simple laboratory or simulation set-ups to facilitate and enhance greater 
understanding of mutual coupling and its effects. 
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I. INTRODUCTION 


A. BACKGROUND 

Antenna arrays are widely employed in both commercial and military 
applications. Consequently, there is much research devoted to enhancing the perfonnance 
of the various array configurations used. In particular, mutual coupling between the 
antenna elements in an antenna array is a potential source of performance degradation. 

Depending on the application, errors due to mutual coupling can be significant. 
Mutual coupling variations between the elements are a source of amplitude and phase 
errors. The presence of mutual coupling distorts phase vectors of radiation sources [Ref. 
1]. This can cause severe degradation of the tracking perfonnance in radar as well as 
increasing the bit enor rate in communication antennas, if it is not properly compensated. 
In effect, this phenomenon introduces a noise floor that precludes synthesis of high- 
quality patterns with very low side lobes or detenninistic pattern nulls. Other possible 
effects appear in signal processing arrays, such as adaptive systems, which can be 
extremely sensitive to small errors due to the non-linear processing involved [Ref. 2], 

For example, in the military, with the unmanned aerial vehicle (UAV) fast 
becoming a popular platform for a broad spectrum of battlefield applications, there is 
now an even greater impetus to better understand and counter the effects of mutual 
coupling. With an increasing payload to handle broader applications, it has become a 
major design challenge to accommodate the large number of electronic warfare and 
communications suites as well as the various other onboard processors on the UAV such 
that the individual components can function effectively on its own, as well as a collective 
whole. For this end-state to be possible, the mutual couplings that result from the 
interactions between the various components must be considered in the design process. 

Mutual coupling is most significant between neighboring elements. It is often 
taken that the radiation pattern of an array of identical antenna elements is the product of 
an element factor and an array factor, on the assumption that all the elements have equal 
radiation patterns [Ref. 3]. For a practical array, this is not entirely true. For example, in a 
linear array consisting of a large numbers of elements, most of the embedded elements 


1 



experience approximately the same electromagnetic environment, thus exhibiting similar 
radiation characteristics. However, a few elements at both ends of the array encounter 
edge effects. This phenomenon becomes more pronounced in smaller arrays, where 
mutual coupling causes each element to see a different environment and consequently has 
a different radiation pattern from its neighboring elements. 

The demand, hence versatility, required of small phased arrays due to the lower 
cost necessitates that the effects of mutual coupling be understood and that compensation 
techniques be developed in the applications of such arrays. 


B. OBJECTIVE OF STUDY 

The vast majority of previous research has dealt with mutual coupling in an 
infinite array environment [Ref. 4-6]. This thesis presents a fundamental quantitative 
analysis of mutual coupling associated with a single dipole and a three-element array in 
their respective simulated operational environments. The mutual coupling effects on the 
impedance of a radiating elemental dipole and on the radiation patterns of the dipole 
array for both transmitting and receiving modes will be investigated. 

The examination of a small array provides a useful environment in which to 
develop, optimize and evaluate the radiating elements. In essence, this study provides a 
better understanding of the effects of mutual coupling and facilitates the design of 
corresponding compensation techniques in practical antenna arrays. 


C. APPROACH 

Commercially available software packages employing such techniques as the 
method of moments (MM) and the finite element method (FEM) are used to study the 
mutual coupling effects of the different dipole array geometries in simulated array and 
ground plane environments. 

Both Patch and NEC-Win Professional are first used to investigate the mutual 
coupling effect of a single dipole that is centrally excited by a voltage. Both tools use the 
MM technique. 
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Next, the mutual coupling effect of a three-element array is studied using Patch. 
The antenna topology under investigation is a linear array of horizontally oriented 
dipoles, each of which is excited by a complex voltage according to the dipole position in 
the array. The scenarios for transmit and receive modes as well as with tenninal loading 
are also investigated. For each of the geometries studied, the impedance and pattern 
results obtained from simulations are compared with the corresponding theoretical 
results. 

D. OVERVIEW 

Chapter II gives a theoretical analysis of the impedance characteristics of a single 
dipole as a radiating element in free space as well as within the confine of a perfect 
ground plane, and of a planar three-element array. This is followed by a discussion of the 
MM as a numerical analysis technique to predict the impedances between elements of an 
antenna array due to mutual coupling. The principle of pattern multiplication in antenna 
radiation and antenna directivity will also be covered. 

Chapter III describes the simulations conducted using the software tools and 
presents the results. The dipole and array under test are subject to different operating 
environments such as geometry and various load conditions. The results from the 
different tools are compared and discussed in relation to mutual coupling and its effects. 

Chapter IV discusses possible compensation techniques that can be explored to 
help reduce the effects of mutual coupling in small dipole arrays. 

Finally, Chapter V summarizes the research and suggests further work with regard 
to compensating for the mutual coupling effects in small dipole arrays. 
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II. THEORY 


A. DIPOLES 

The properties of any array depend on the characteristics of the individual 
radiating elements. The dipole in its various forms constitutes one of the most widely 
used radiating elements in arrays. This section begins with a discussion of a dipole in free 
space, followed by an analysis of its characteristics in the presence of a perfect electric 
conducting (PEC) ground plane to illustrate the effect of mutual coupling. 

1. Single Dipole in Free Space 

The fundamental building block of a practical dipole is the ideal Hertzian dipole. 
The ideal dipole is an infinitesimal element with a current of unifonn magnitude and 
phase. Therefore, the radiation field from a longer practical length dipole is essentially 
the vector sum (integration) of the contributions from all ideal dipoles weighted by the 
current distribution. 

The dipole field pattern in an array is called the element pattern, EP , which is 
exactly sin 0 for the Hertzian dipole. On a finite dipole, the current amplitude and phase 
are no longer constant. For a short dipole where its length L«/ 1/2, the current is in 
phase but its amplitude is approximately triangular. For a half-wave dipole, the current is 
still in phase and its amplitude can be approximated by a sinusoid. For these dipoles, the 
radiation will be strongest in the direction normal to the dipole and weakest along the 
axis of the dipole. As the dipole length increases, the radiation pattern sharpens, until at a 
length of A/2 , the half-power beam width narrows from 90° to 78°. However, for longer 
dipoles, the current on some sections will be out of phase with others, leading to partial or 
total cancellation in the far field [Ref. 3] 

The infinitely thin A/2 dipole in free space has a center-fed radiation resistance 
of 73.1 Q. At dipole length L = A/2 , the impedance is slightly inductive and the total 
radiation impedance is 73.1 + J42.5 £1 . This inductive impedance drops rapidly to zero as 
the dipole is foreshortened. For L = 0.485/1, the dipole is resonant and the input 
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impedance is 71.8 D.. Below resonance, the reactive component of the dipole impedance 
is capacitive. In the range below A/2, the dipole radiation resistance decreases 
monotonically with length and is almost independent of diameter, whereas the reactive 
component depends heavily on diameter. The thinner antennas are more capacitive for a 
given length [Ref. 1]. 

In practice, a dipole has a finite length-to-diameter ratio. Additionally, the 
environment also affects its impedance. Both the radiation pattern and impedance are 
influenced by the presence of nearby objects. The most commonly encountered object is 
the ground plane, which is intentionally introduced to increase the dipole’s directivity. 

2. Single Dipole Above a Perfect Ground Plane 

The dipole impedance will change if the antenna environment is varied. As the 
impedance changes, the current on the dipole may become redistributed and thus alter the 
dipole’s radiation pattern. 

The presence of a ground plane affects the dipole’s impedance. Consider a single 
dipole above an infinite ground plane that is also a perfect electric conductor. The image 
theory can be used to derive the relationships between the excitation voltage, impedance 
and current of the dipole. The image-method equivalent of a current-carrying dipole 
above a PEC ground is shown in Figure 1. 
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a. Horizontal: 



T I 


Figure 1. Image equivalent of a single dipole above ground. 

In each case, the voltage across the dipole is governed by the following 
expression: 

V x = I\Z n +I 2 Z l2 (III) 

where/j and /, are the source and image currents respectively, and Z n and Z 12 are the 
self and mutual impedances respectively. 

Applying the image current relations to Equation (II. 1), the voltage and active 
impedance across each dipole become as follows: 

For the horizontal dipole, 


V — T 7 -T 7 
v \ 1 \^\\ 1 1^12 

(11.2) 

„N 

II 

| 

II 

JN 

1 

JN 

K> 

(II.3) 

For the vertical dipole, 


v = I 7 +17 

(II.4) 
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(II.5) 


V 

7 — — = 7 +7 

Z.J Zs n -rzs l2 


Equations (II.3) and (II.5) give the impedance of a horizontal and vertical dipole 
over a ground plane respectively. In each case, the impedance is a function of the mutual 
impedance, Z 12 , which is the result of mutual coupling between the dipole and the 
ground. By reciprocity, the mutual impedance of the image dipole is equal to that of the 
source dipole itself, i.e. Z 21 = Z 12 . 


The image method in effect produces a pair of parallel dipole antennas, side-by- 
side in Figure la and collinear in Figure lb respectively. Note that in order to apply the 
image method, the ground plane must be infinite. However, the method is often applied 
to finite ground plane if the edges extend sufficiently beyond the array of dipoles. 

For the parallel side-by-side dipole antenna, the mutual resistance and mutual 
reactance are [Ref. 7] 


R n = 3012cos(/?/?)-cos 


Ih +L +L 


■cos 




Ir + L +L 


(II.6) 


X l2 = — 3oj2sin(/?A)-sin /3[^! 


h +L +L 


-sin 


fi^lr+L 2 +l) } (II.7) 


where (3 = ■ 


For the parallel collinear dipole antenna, the mutual resistance and mutual 
reactance are [Ref. 7] 


R n — -15 cos (j3h) 


-2 cos (2jBh) + cos (2/3 (h - L)) + cos (2/5 (h + Z)) - In 


r h 2 -L 2 ^ 
V h 2 j 


+15sin(/?/2)[^2sin(2/?/i)-sin(2/?(/?-Z))-sin(2/?(/?-l-Z)) (II. 8) 


X n = 15sin(/?/2) 


2 cos (2j3h) - cos (2y? (/? - L )) - cos (2y9 (/? + L )) - In 


r h 2 -L^ 

V h 2 j 


-15cos(/?/2)[2sin(2/?/2)-sin(2/?(/2-Z))-sin(2/?(/2 + Z)) (H-9) 
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where in both cases, 


Z,2 — ^12 + 7^12 - % 2 \ — ^21 + 7^21 (II. 10) 

In the presence of other dipoles, such as in a dipole array, mutual coupling will be 
more complex. Consequently, an array designer must consider the theoretical and 
empirical aspects of matching the dipole in an array environment. 


B. DIPOLE ARRAYS 

Dipole arrays are rarely designed to radiate into full space - the dipole elements 
are usually backed by reflectors, such as earth or wire grids. As the more central elements 
in an array are negligibly affected by the ground-plane edges, they can be viewed as 
being backed by an infinite ground plane. However, the few elements at both ends are 
significantly affected by mutual coupling, even if the ground plane is infinite. 

To better understand this effect, consider a three-element dipole array parallel to 
and above an infinite ground plane. The image theory is again used to derive the 
relationships between the excitation voltage, impedance and current of the individual 
dipole. The image-method equivalent is shown in Figure 2. 



Figure 2. Image equivalent of a 3-dipole array above ground. 

By the image method, in general, 

h=~h, h=~h, h=~h ( 11 . 11 ) 

By symmetry, at broadside scanning, 

h=h> h=h=~h (H-12) 
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The currents in all the source dipoles as well as their image dipoles contribute to 
the voltage across each dipole element. For dipole 1, the voltage across it is given by the 
following expression: 

V \ =I i Z n + / 2 Z| 2 + /3Z13 + 1 4Z 14 + f 5 Z| 5 + h Z \6 


■ /jZj, + I 2 Z l2 + 7 3 Z 13 • 


h Z 14 


-/ 2 Zi5 + 7 6 Z 16 


— 7i(Zn Z 14 ) + 7 2 (Z 12 Z 15 ) + 7 3 (Z 13 Z 16 ) 

= 7 1 z; 1 +7 2 z; 2 +7 3 z; 3 ( 11.13) 

where Z' is the mutual impedance of dipole 1 with dipole /' and its image. 

By the same account, for dipoles 2 and 3, 

V 2 = I\Z'ix + 7 2 Z 22 + 7 3 Z 23 (11.14) 

F 3 =7 1 z; 1 +7 2 z; 2 +7 3 z 3 ' 3 (11.15) 

The corresponding active impedances for the dipoles are derived as follows: 

F II 

Z = 7 J- = Z' + 2 ^ 7 ' + —Z' 

^1 j ^11 ~ j ^12 T j ^13 

■'l -*1 A 

z 2 =f=fz:,+x+fz: 3 

Z, = Y = yZ' s , + ^Z' s , + z: s 

^3 ^3 ^3 

Equations (11.13) to (11.15) can be represented in matrix form as 

[F] = [Z'][7] (11.19) 

For more complex geometries, such a simple analysis will not be possible. In 
place, numerical methods are used; most commonly, the Method of Moments (MM), and 
Finite Difference Method. The MM is the primary method used in the analysis and is 
discussed in the next section. 


(11.16) 

(11.17) 

(11.18) 
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C. METHOD OF MOMENTS 

The MM is a numerical technique that is applicable to arbitrary bodies. 
Essentially, it reduces the E-field integral equation (EFIE) into a matrix problem, the size 
of which is related directly to the electrical size of the antenna in terms of wavelength. 
Large bodies result in large matrices and, therefore, require large computers with fast 
processing units. The use of MM in the analysis and design of arrays of wire antennas (or 
scatterers) has significant advantages over the more classical methods used in treating 
arrays in that mutual coupling between array elements is taken completely into account. 
Furthermore, realistic assumptions can be made of the current distributions on the wires 
and the type of wire element array problem that can be considered is rather general, given 
that the array elements can be excited or loaded at any points. The use of MM to predict 
the impedances due to mutual coupling is discussed in the ensuing paragraphs [Ref. 8]. 



Figure 3. Derivation of EFIE. (After: [Ref. 8]) 

The first step in applying the MM is to represent the current on an arbitrary 
perfect electric conductor (PEC) by a series of unknown expansion coefficients I n : 

= ( 11 - 20 ) 

n=l 

where J n are the basis functions, and both I n and J n can be complex. The selection of 

the basis functions is an important consideration in the MM solution. Basis functions 
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should be mathematically convenient and be consistent with the behavior of the current. 
The series representation for the current density is then inserted into the EFIE. 

The EFIE for a PEC is [Ref. 8] 


£.(OL = £/„JT W/„(r)G(F,r)-^-[v'-y.(?)lv'G(r,F)l ds 

»=1 a> £L J J tan 


( 11 . 21 ) 


where co is the angular frequency (assuming a time dependence e jM ), fl the 
penneability, £ the permittivity, r a position vector to an observation point on the 
surface of the PEC and G(r,r') the free space Green’s function denoted by 


e ~ ] P R 

G {r,r) = — (11.22) 

4 kR 


Next, define a set of weighting functions using Galerkin’s method [Ref 8]. The 
weighting functions are chosen as follows: 


W, =J 


(11.23) 


where m = 1,2, ... N. 


Both sides of Equation (11.21) are multiplied by each of the weighting functions 
and integrated over its domain, deriving N equations of the form: 

JL W ,n O 7 )' E t (?) | tan ds = jj Sm W m (r) ■ j j-\/„ x JJ jO)flJ n (r')G(r, r) 


CO£ 


V' ■ J n (C)J V'G(r, r) ds'} ds 

J tan 


(11.24) 


Now, define the impedance element: 


ff ds \ 

1_ 

^3 

JJSm J 

J Sn 


jafiw m (r) ■ J n (r)-My'- J n (?'m ■ w m (f )) 

COE 


G(r,r') (11.25) 


where the surface divergence theorem has been applied: 

JJX- v 'g*=£g(v^ 


(11.26) 
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Let the right-hand side of Equation (11.24) be defined as 

K=\\ Sm Kir)-E,(f)ds (11.27) 

Equation (11.24) can now be written as 

V =Y/Z (11.28) 

m / j n mn v ' 

n =1 

or in matrix fonn: 

K] = [ZJW (H.29) 

The matrix V is the excitation vector and Z is the impedance matrix. The vector/ 

m mn f n 

contains the unknown expansion coefficients, which are determined by solving: 

= (11.30) 

Once these expansion coefficients are determined, the series representation for the 

current, Equation (11.20), is defined and can therefore be used in the radiation integral to 
obtain the electric field at the observation point. 

The MM is a powerful tool that provides rigorous solution for the induced current 
density on a body. Various proofs and theorems assure that the MM solution is rigorous 
under realistic conditions. In general, if the MM is properly applied, the series 
representation for the current will converge to the actual current as the number of basis 
functions is increased. 

D. PATTERN MULTIPLICATON AND DIRECTIVITY 

By the principle of pattern multiplication, the field pattern of an array consisting 
of similar elements is the product of the pattern of one of the elements (element pattern, 
EP ) and the pattern of an array of isotropic point sources with the same locations, 
relative amplitudes, and phases as the original array (array factor, AF ) [Ref. 3]. 
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Figure 4. Single dipole parallel to and above perfect ground. 


Consider a dipole of length L parallel to and of height h above a perfect ground 
plane, as shown in Figure 4. The normalized element pattern of the single dipole in free 
space is [Ref. 3] 


EP ■ 


i i/2 

- J e i/}z '™ e dz -- 
L i/2 


: sin# 


sin (/?T/2)cos#] 
(J3L/ 2)cos# 


(11.31) 


From Figure 4, the far field phase difference between the dipole and its image is 
2 jsm(fih cos a ), where a is the scan angle of an incoming (or outgoing) plane wave 
relative to the nonnal (y -axis) to the ground plane. The dipole and its image in effect 
form an array of two dipoles, with the array factor 

AF = 2 j sin (/3h cos a) 

= 2 j sin(/?/z sin </>) (11.32) 


where a = 90° -(/>,(/> being the angle measured from the x -axis in the x-y plane. 

By the principle of pattern multiplication, the overall field pattern F for the 
configuration in Figure 4 is therefore 

sinr(/?L/2)cos#l 

F = EPxAF = sin#—--- -2 jsm(/3hsm </>) (11.33) 

[pL/ 2)cos# 

Next, consider an equally spaced linear array of N identical dipoles, each dipole 

with a current of unit amplitude and uniform phase, at a height h above a perfect ground 

plane as shown in Figure 5. For this configuration, the element pattern consists of that for 

a dipole above the ground plane as defined by Equation (11.33), with the array factor 
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consisting of N isotropic point sources with the same locations, relative amplitudes, and 
phases as the original array. 



Figure 5. Equally spaced linear array of N identical dipoles above perfect ground. 


From Figure 5, the phase difference between any two adjacent dipoles is 

y/ = /3dsina = /3dcos0 (11.34) 


The array factor is given by [Ref. 3] 

A T7 = V Jw = J( N ~ AvI 2 si n { N V/2) 
h sin(^/2) 


When normalized, AF becomes 


AF = 


sin(iV^/2) 

Nsin(i/f/2) 


(11.35) 


(11.36) 


Similarly, by the principle of pattern multiplication, the complete normalized field pattern 
F of the antenna array is then given by 


F = sin#- 


sin 


[W 2) 


COS 


0~\ n smiNy/12) 

- -—-^ 2 j sin (J3h sin </)) - v Yl { 

(y?E/2)cos6> A^sin(^/2) 


(11.37) 


Due to the ground plane, this field is radiated only over a half sphere above the 
ground plane. As a measure of the radiation intensity in a certain direction relative to an 
isotropic antenna, the directivity of this array can be derived. By definition, the directivity 
is given by 


D{0,4>) 


An 

^A 


(11.38) 
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where Q A is the beam solid angle defined by 


n A (W) = II sindddd(/> (11.39) 

0 0 

To summarize, in this section, the radiation pattern of an array of dipoles over an 
infinite ground plane is derived assuming identical element patterns. In practice, due to 
mutual coupling, the total array pattern, hence directivity, deviates from the ideal case. In 
addition, the mutual coupling effects depend on the frequency and scan direction. The 
simulations in the next chapter show the effects of mutual coupling. 
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III. SIMULATIONS AND RESULTS 


A. SIMULATIONS 

1. Simulation Parameters 

The frequency used in the simulations is 300 MHz, which results in a one-meter 
wavelength. The basic dipole geometry consists of a quadrilateral of length L = 0.45A 
and width w = 0.005/1. This length is chosen so that the current on the entire dipole is in 
phase, like in most practical applications. The length-to-width ratio is 90, which closely 
approximates an infinitely thin dipole. 

2. General Procedures 

The antenna analysis software packages, Patch and NEC-Win Professional, are 
used to determine the active impedance of a centrally voltage-excited dipole in free space 
and its mutual impedance with a perfect ground. The dipole is placed parallel as well as 
perpendicular to the ground plane as in Figure 1 and its distance above the ground plane 
is varied. The results obtained are compared with those from Equations (II.6) to (II.9). 

Next, the effect of mutual coupling on the phase and directivity of a dipole array 
over a perfect ground plane is investigated using Patch. The array topology comprises a 
linear array of three horizontally oriented dipoles, each of which is excited by a complex 
voltage according to the dipole position in the array, as shown in Figure 6. The scenarios 
for transmit and receive modes as well as with terminal loading are also investigated. 



d sin a 


•v< 


.jjjj* 


/ dipole 1 dipole 2 dipole 3 

Figure 6. Three-dipole array above a ground plane. 
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3. Cylindrical Equivalence of a Quadrilateral 

The quadrilateral dipole model is used in Patch. In NEC-Win Professional, a 
cylindrical dipole model is used instead. The equivalent cylindrical dipole has a length 
L = 0.45 X and a radius r = 0.00125/1, where its radius is related to the width of the strip 
by r = 0.25w [Ref. 9]. 


4. Convergence in Patch 

Patch uses triangular sub-domains in the MM formulation. The patch geometry 


for the three-dipole array over a ground plane is depicted in Figure 7. 



For convergence, the segment length of each triangular patch should be less than 
0.1 X. The smaller the segments, the more accurate the result will be. However, there will 
be a corresponding increase in computation time. 

B. RESULTS 

1. Impedance of a Dipole in Free Space 

Using Patch, the current density at the center of the dipole is found to be 
J = 1.7876+ /'l.7136 A/m, giving a total current I = Jxw = 8.938 + y'8.586 inA. With a 
unit voltage excitation, the self-impedance of the dipole is calculated to be 
Z n = 58.3-./55.89 £2. 

Using NEC-Win Professional, the corresponding impedance obtained for eleven 
segments is 61.803- /43.398 U . Slight variations were then made to both the length and 
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radius of the cylindrical dipole. It was found that the closest agreement with Patch occurs 
when the cylinder length is shortened from 0.45A to 0.44/1. At this length, the 
impedance becomes 58.24-y'58.89 £2 . Although Patch and NEC-Win Professional use 
the same numerical technique for antenna analysis, differences inherent in the coding and 
assumptions are expected. This length of 0.44/1 for the cylindrical dipole in NEC-Win 
Professional will be used in subsequent simulations. 

As expected, the impedance obtained is capacitive, since the dipole length is less 
than 0.485/1, the approximate length at which resonance occurs. This result also agrees 
with the general behavior of input impedance based on actual measurements carried out 
by Brown and Woodward [Ref. 10]. 

2. Mutual Impedance of a Dipole Above Perfect Ground Plane 

For each height above the ground plane, the impedance of the dipole Z, is 

obtained as in Section III.B.l. To derive the mutual impedance Z 12 , its self-impedance 
Z n is subtracted from its active impedance, Z x , using Equation (II.5). The plots of Z 12 
against the height above the ground plane for a horizontal and vertical dipole are given in 
Figures 8 and 9 respectively. 



Figure 8. Mutual impedance of a horizontal dipole above ground plane. 
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Black: Theory Blue: Patch Red: NEC-Win Pro 



Figure 9. Mutual impedance of a vertical dipole above ground plane. 


The mutual coupling between the dipole and the ground plane introduces Z 12 . 
The plots in Figures 8 and 9 show that as the dipole is moved further away from the 
ground plane, both the resistive and reactive components of Z 12 approach zero. This is 
expected because the coupling influence of the plane diminishes with increasing height, 
thus approximating a free space condition. 

The mutual impedance obtained from Patch closely agrees with that from NEC- 
Win Professional. Both sets of simulation results are also in line with the corresponding 
theoretical values given by Equations (II.6) to (II.9). The results further establish the 
validity of MM as an analytical tool to determine the mutual coupling in an array. Since 
Patch and NEC-Win Pro give similar results, only Patch is used in subsequent analyses. 

3. Effect of Mutual Coupling on 3-Dipole Array Above Perfect Ground 
Plane 

The effect of mutual coupling on the phase, directivity and active element pattern 
of the 3-dipole array is investigated for dipole separation of 0.5A and 0.375/1, with the 
array at a height 0.25A above the infinite PEC ground plane. In the transmit scenario, 
each dipole is centrally excited by a complex voltage of unit amplitude according to the 
dipole position in the array when the array is scanned. In the receive scenario, an 
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incoming plane wave of unit amplitude is used. Thus, in the ideal situation, the phase 
should be linear and the amplitudes equal. 


a. Phase 

The effect of mutual coupling on the terminal phases of the array dipoles 
as a function of scan angle a is investigated for both transmit and receive modes. 


dipole separation = 0.5 wavelength 



Figure 10. Phase difference on transmit as a function of scan angle a for d = 0.5/1. 


dipole separation = 0.375 wavelength 



Scan angle (degrees) 


Figure 11. Phase difference on transmit as a function of scan angle a for d = 0.375/1. 
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In Figures 10 and 11, the blue curve is obtained by subtracting the 
terminal phase of dipole 1 from that of dipole 2 and the green curve is obtained by 
subtracting the terminal phase of dipole 2 from dipole 3. The red curve is the theoretical 
value for both the blue and green curves, as governed by Equation (11.34). 

In theory, the tenninal phase of dipole 3 leads that of dipole 2, which in 
turn leads that of dipole 1. This means that for scan angles other than broadside, both the 
blue and green curves should be greater than zero. The results in Figures 10 and 11 show 
that for small scan angles (less than 8° and 12° for the respective arrays), there is a phase 
lag for the green curve instead. At greater scan angles, the phases are also found to 
deviate from Equation (11.34). Both sets of results exhibit similar behavior. 

The phase deviations between adjacent dipoles relative to the theoretical 
value are derived based on Figures 10 and 11 and plotted in Figures 12 and 13 
respectively. 


dipole separation = 0.5 wavelength 



Figure 12. Phase deviation on transmit as a function of scan angle a for d = 0.5/1. 
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dipole separation = 0.375 wavelength 



Figure 13. Phase deviation on transmit as a function of scan angle a for d = 0.375/1. 


From Figures 12 and 13, it is observed that for both arrays, the phase 
deviations are greatest at broadside. For Figure 12, the deviation is minimum at scan 
angle of about 20° whereas for Figure 13, this occurs at about 16°. Beyond these angles, 
the deviations fluctuate but appear to smoothen out towards the array endfire. 


dipole separation = 0.5 wavelength straight:transmit dotted: receive 



Figure 14. Comparison of phase difference between transmit and receive. 
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The phase difference between adjacent dipoles as a function of scan angle 
for d = 0.5A on receiving incoming plane waves is also investigated and compared with 
the corresponding values on transmit. The results are plotted in Figure 14. The two sets of 
results agree closely. This observation implies that the reciprocity theorem holds even in 
the presence of mutual coupling. The theorem therefore remains a very useful design 
tool. 

As a result of mutual coupling, mutual impedances between the dipoles 
have been introduced upon excitation of the dipoles. This in turn distorts the phases at 
each dipole terminal, accounting for the deviations observed in Figures 10 to 13. 
Concomitantly, the directivity of the array is expected to vary with scan angle. The effect 
on directivity will be investigated in the next section. 

b. Directivity 

The effect of mutual coupling on array directivity in the transmit mode as 
a function of scan angle a is investigated. 


dipole separation = 0.5 wavelength 



Figure 15. Array directivity on transmit for d = 0.5/1. 
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dipole separation = 0.375 wavelength 



Figure 16. Array directivity on transmit for d = 0.375/1. 


From Figure 15, for d = 0.5 A , the realized array directivity is less than the 
theoretical values obtained from Equations (11.38) and (11.39) at all scan angles. The 
difference is minimum at scan angle of 15° and maximum at scan angle of 50°. At 
broadside, the difference is about 0.7 dB. From Figure 16, for d = 0.375/1, for most scan 
angles, the realized directivity is also less than the theoretical values. For this array, the 
difference at broadside is about 0.25 dB 

The assumption made in arriving at Equation (11.39) is that the element 
terminal currents are proportional to their excitations, the current distributions on the 
elements in the array are identical and pattern multiplication is valid. However, in a real 
array, the elements interact with each other and alter the currents, thus impedance, from 
that which would exist if the elements were isolated. This interaction, or mutual coupling, 
changes the magnitude, phase and distribution of current on each element and in turn is 
manifested in the loss of directivity of the arrays, hence the results obtained in Figures 15 
and 16. 
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c. Active Element Pattern 

The active element pattern of three-dipole array is obtained by centrally 
exciting a single element in the array with a unit voltage, with the other two elements 
terminated in loads of 50 £2. Three cases are investigated: (1) all the dipoles are centrally 
excited, (2) the center dipole, i.e. dipole 2, is centrally excited and (3) a side dipole, 
dipole 3, is centrally excited. The pattern plots for dipole separation of d = 0.5/1, 
d = 0.375/1, d = 0.25 A are given in Figures 17 to 19 respectively. 


theta = 90 degrees dipole separation = 0.5 wavelength 



Figure 17. Active element patterns for d = 0.5 A . 


theta = 90 degrees dipole separation = 0.375 wavelength 



Figure 18. Active element patterns d = 0.375/1. 
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theta = 90 degrees dipole separation = 0.25 wavelength 



Figure 19. Active element patterns for d = 0.25 X. 


In Figures 17 to 19, the red and green curves represent the active element 
patterns of the array as defined by Case (2) and Case (3) respectively. These patterns are 
different from what is expected of an isolated radiating dipole above a perfect ground, as 
defined by Equation (11.33). This is due to mutual coupling between the radiating dipole 
and its adjacent dipoles, which causes currents to be induced on the adjacent dipoles. 
Consequently, the adjacent dipoles also radiate power, albeit at a lower level compared to 
the power radiated from the excited dipole. In Case (2), because the excited dipole is 
centrally located, the induced currents, hence impedances, on dipole 1 and dipole 3 are 
equal (Equations (11.12), (11.16) and (11.18)). Consequently, the green curve is 
symmetrical about broadside, as opposed to the skewed effect observed in the red curve. 
In addition, the red curve shows gain reductions at scan angles of about 35°, 25° and 15° 
for d = 0.5A, d = 0.375/1 and <7 = 0.25/1 respectively. At these corresponding scan 
angles, Figures 15 and 16 on directivity show comparable drops in gain when scanned to 
these angles. This observation agrees with the theoretical analysis made by Pozar [Ref. 
12 ]. 

The blue curve represents the gain contributed by all the three radiating 
dipoles in the array and the corresponding mutual coupling effects. As a result, its gain at 
broadside is higher than those of the other two curves. In theory, at broadside, gain is 
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proportional to the array dimension, thus accounting for the observations made where the 
corresponding gain is highest for d = 0.5 A and lowest for d = 0.25A . However, side 
lobes exist for the arrays with d = 0.5/1 and d = 0.375/1. Ideally, nulls (zero gain) would 
exist between the main lobe and sidelobes, but mutual coupling causes a filling in of the 
nulls. 


C. DISCUSSIONS OF RESULTS 

The effect of mutual coupling diminishes as the separation between dipoles 
increases, i.e. it is most significant between neighboring elements, and this effect is more 
pronounced when the dipoles are placed side by side than when they are collinear. These 
observations are also intuitively true. In the fonner, a large separation approximates a 
free space condition, and in the latter, the effective surface area for interaction is much 
reduced for collinear dipoles as compared to side-by-side dipoles. The driving point 
impedance, or active impedance, is found to be a function of terminal currents and mutual 
impedance. 

In a scanned array, the phases of the tenninal currents are also varied under the 
influence of mutual coupling, and the errors are found to be greatest at broadside and 
appear to smoothen out at the array endfire. This effect of mutual coupling in an array is 
manifested by the far field pattern, or more specifically in the simulations, the directivity 
of the array. As opposed to the classical theoretical approximation obtained by 
considering each element in isolation, the actual realized directivity include variations in 
the excitation currents as well as the patterns of each element acting under the influence 
of all coupling effects. Mutual coupling changes the current magnitude, phase and 
distribution on the dipole, thereby ‘distorting’ the radiation vectors, giving rise to peaks 
and nulls in directions that are different from theoretical approximations. 

The radiation properties of an array can be evaluated using the approach of active 
element pattern, which is obtained by exciting only an element while match-loading all 
other elements in the array. The obtained active element pattern is then due to the direct 
radiation from the excited element combined with the fields re-radiated from the other 
elements, which in turn receive their power through spatial coupling with the excited 

element. This coupling is a function of element characteristic and the array geometry. 
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The active element pattern depends on the position of the fed element in the array, so that 
the edge elements will have different active element patterns from those of the elements 
near the center of the array. If the array is large, however, most of the elements will see a 
uniform neighboring elements environment and the pattern can be approximated as equal 
for all elements in the array. To represent the possibility of gain variations, the active 
element levels are relative to a reference element located at the center of the array. 

The importance of the results on active element pattern comes from the fact that 
measurements of the active element patterns can give more directly the desired scan 
characteristics of an array under test and are therefore useful in predicting the scan 
performance of large phased arrays, as illustrated by the close relation between nulls in 
the active element patterns and losses in directivity. 
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IV. COMPENSATION TECHNIQUES 


A. INTRODUCTION 

Arrays are now widely used in electronic warfare and communications 
applications. Mutual coupling between elements can impact the array performance. The 
mutual coupling tenns for elements near the edges of an array are significantly different 
from those for elements in the interior of a large array. This edge effect is particularly 
significant for small arrays. 

The simulations have shown that the active impedance of each dipole element in 
an array varies as a function of the beam scan angle due to mutual coupling. In general, 
this results in a loss in array directivity, hence gain. The reduction in gain degrades the 
performance of the array, thereby limiting its usefulness. For example, for small direction 
finding arrays, mutual coupling leads to large phase errors that degrade the angle of 
arrival estimates, while for large arrays that have ultra-low sidelobe amplitude 
distributions, mutual coupling can affect the sidelobe levels, even though the majority of 
elements are away from edges. 

Much research work has gone into the compensating techniques for various array 
configurations. Each technique has its own advantages and disadvantages. In the context 
of the simulations done in this work, some possible compensation techniques so that the 
array can be steered over a wide scan angle without significantly affecting its 
performance are suggested and discussed in the ensuing section. 

B. COMPENSATION TECHNIQUES 

1. Modification Within the Feed Line 

The goal of this technique is to reduce the input impedance variation with scan 
angle. Theoretically, this can be accomplished by adding loads to the transmitters or 
receivers via series impedance or terminated circulators, as shown in Figure 20. The 
compensating impedance Z c (i = 1,2,3) adds to the corresponding impedance Z ; in 

Equations (11.16) to (11.18). Thus, variations in the impedances from element to element 
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can be corrected. However, the same Z c must be used for all angles and therefore only an 

average correction can be introduced. The circulators do not correct for mutual coupling 
per se, but isolate the transmitter from reflections, which change as a function of angle 
due to mutual coupling, from the elements. For this technique, the concomitant effects 
such as the intrinsic loss due to the modification and variation of realized gain and 
bandwidth with scan angle should be investigated in tandem. 


a. Series impedance b. Terminated circulators 



Figure 20. Techniques to impedance-match elements. (After: [Ref. 1]) 


2. Modification in External Environment of Array 

This is based on the concept originally suggested and investigated by Edelberg 
and Oliner [Ref. 1]. One approach to mutual coupling mitigation is a simple add-on of 
dummy elements, which are loaded elements at the edges of the array, as shown in Figure 
21. In a sense, the dummy elements “soften” the edge. 


frequency = 7.65 GHz 
/dummy 


v ^ 
A 71 


dummy 7 


•v</ 


L 0-452 ^/o.506^/o.5Q6^/ 0,45 kj /\ 

WBm rv ' 


dipole 1 dipole 2 dipole 3 


/ground 

Figure 21. Three-dipole array with a dummy element on each side. (After: [Ref. 11]) 


Based on this configuration, preliminary results have been obtained using Patch 


and are given in Figures 22 and 23 [Ref. 11]. The data was generated by taking the 
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transmit active element patterns and using them as receive element patterns when all 
elements are combined on receive. Ideally, the phase across the array should be linear. 
The linear phase is subtracted out and the residual errors are plotted in the figures. Using 
dummy elements with a 50 £2 load, the maximum phase error has been reduced by about 
half. However, this load is not necessarily the optimum. For this technique, the number of 
dummies, their configurations and the optimum load conditions to closely approximate an 
infinite array environment need to be determined. 


error in phase difference 



Figure 22. Phase error for the array without dummy elements. (From: [Ref. 11]) 


error in phase difference 



Figure 23. Reduced phase error for the array with dummy elements. (From: [Ref. 11]) 
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3. Connecting Circuits 

In this technique, matching circuits that vary with scan angle are used, so as to 
impedance-match phased arrays over wide scan angles. Such a conceptual network is 
shown in Figure 24 for collinear dipoles. Multi-port elements are used, with one of the 
feed points serving as a source for coupling corrections. A candidate element design that 
has shown promise for this application is a hybrid fed printed circuit dipole. It allows the 
adjustment of the free space coupling signals to be cancelled by the compensation 
network [Ref. 11]. The advantage of this approach is that it provides compensation that 
varies with angle. Hence, it has the potential to effectively correct for mutual coupling, 
which also varies with angle. 


multi-feed point elements 



Figure 24. Compensation network integrated into a multi-feed array. (After: [Ref. 11]) 

The matching circuit must be connected in such a way as to preserve the 
symmetry of the array. In a planar array with identical elements, this simply means that 
the connecting circuits should be added in sets of parallel rows, each having identical 
circuits. Likewise, this technique attempts to achieve an infinite array environment. There 
are various forms that the connecting circuits may take. Capacitors, inductors, 
transmission lines or even more complex active networks can be employed, the choice 
being dependent on practical considerations, which vary from one antenna to another. To 
adopt this technique, the analysis, design and simulation of a mutual coupling 
compensation network and its integration with the beam-forming array need to be 
conducted. 

With the requirements for operational arrays becoming more stringent and 

demanding, it is inevitable that for greater operational effectiveness, the development of 

compensation techniques to counter the effects of mutual coupling has also become an 
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integral design component. In practice, it is extremely difficult to completely control 
mutual coupling. Additionally, complexity and cost are also major design considerations. 
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V. CONCLUSIONS 


This thesis work has established the effects of mutual coupling between a 
radiating dipole element and its environment, and between radiating dipole elements 
within a 3-dipole array in a simulated operational environment. The simulations have 
shown that the active impedance of any element in an array depends on both the self and 
mutual impedances and that the terminal current phases, hence directivity, of an array are 
affected by mutual coupling such that the array pattern deviates from that defined by 
classical theoretical approach. In an active element pattern approach, the mutual 
couplings in an array are accounted for through the active element, making this approach 
a viable one to study the radiation patterns of practical dipole arrays. Although the results 
have only been obtained for a single dipole and 3-dipole array above a simulated ground 
plane, they are based on the use of MM, which accounts for mutual couplings and is 
generally applicable to any geometry. The results can therefore be extended to arrays of 
different geometries and configurations. The examination of a small array also provides a 
useful environment in which to develop, optimize and evaluate the radiating elements. 

It should be noted that the presence of mutual coupling is not necessarily bad; in 
fact, Hannan has shown that mutual coupling is actually required in closely spaced 
infinite arrays to achieve idealized scanning performance [Ref. 13]. However, the 
presence of mutual coupling variations in a finite antenna array causes errors, the extent 
of which depends on the geometry and configuration of the array. Depending on the 
applications, such errors can be significant and result in severe system performance 
degradation, and must therefore be compensated for accordingly. In light of this, the 
research has also discussed further works with regard to compensating for the mutual 
coupling effects in dipole arrays. Even if these techniques do not result in the desired 
end-state of implementation in operational arrays, they can still be useful in simple 
laboratory or simulation set-ups to facilitate and enhance greater understanding of mutual 
coupling and its effects. 

The next phase of this research should focus on the compensation techniques of 
dummy elements and compensation networks. Dummy elements are appealing because 
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they are simply an add-on solution that does not require redesign of an existing array. The 
optimum number, spacing and load conditions must be determined. 

The compensation network depicted in Figure 24 would have to incorporate some 
complex weights (attenuators and phase shifters). The entire compensation network can 
be modeled using a combination of method of moments and scattering parameters as 
described in [Ref. 14]. This approach includes all the interactions between the elements, 
both via free space (i.e. mutual coupling) and through the feed network. It is conceivable 
that for a small array, a relatively compact solution for the compensation weights might 
be possible. 

In conclusion, this study provides a fundamental understanding of the effects of 
mutual coupling in antenna arrays and facilitates the design of corresponding 
compensation techniques in practical arrays. 
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APPENDIX I. PATCH FILES 


A. INPATCH FILE FOR CENTRALLY EXCITED ARRAY 

A sample inpatch file for the centrally-excited three-dipole array at a scanning 
angle of 15° is as follows: 


(Geometry) 

66 123 

1 0.000000E+00 2.475000E-01 -2.250000E-01 

2 0.000000E+00 2.525000E-01 -2.250000E-01 

3 0.000000E+00 2.475000E-01 -1.800000E-01 


65 0.500000E+00 2.475000E-01 2.250000E-01 

66 0.500000E+00 2.525000E-01 2.250000E-01 

(Edges,faces,vertices) 

1 1 2 

2 1 3 

3 1 4 


122 64 66 

123 65 66 

(Simulation parameters) 

1 

0-10 

1 

v 

3 

21 14 .10000E+01 .00000E+00 

103 58 .68720E+00 -0.7264E+00 

62 36 .68720E+00 0.7264E+00 

0 
0 

.false. 0 0 

1 

90.000000 90.000000 1 0.000000E+00 180.000000 

181 
0 

.true. 

0 

3.000000E+08 

- 1.000000 
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B. INPATCH FILE FOR PLANE WAVES INCIDENT ON ARRAY 

A sample inpatch file for plane waves incident on the three-dipole array at 
different scan angle is as follows: 


(Geometry) 

66 123 

1 0.000000E+00 2.475000E-01 -2.250000E-01 

2 0.000000E+00 2.525000E-01 -2.250000E-01 

3 0.000000E+00 2.475000E-01 -1.800000E-01 


65 0.500000E+00 2.475000E-01 2.250000E-01 

66 0.500000E+00 2.525000E-01 2.250000E-01 

(Edges,faces,vertices) 

1 1 2 

2 1 3 

3 1 4 


122 64 66 

123 65 66 
(Simulation parameters) 

1 

0-10 

180 

P 

90 1 1 0 0 0 

0-10 
P 

90 2 1 0 0 0 

P 

90 178 1 00 0 

0-10 


P 

90 179 1 0 0 0 

0-10 
P 

90 180 1 00 0 

0-10 
0 

0 

.false. 0 0 
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1 

90 90 1 90 90 

1 

0 

.true. 

3 

21 

103 

62 

3.000000E+08 

- 1.000000 
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APPENDIX II. MATLAB PROGRAMS 


A. IMPEDANCE OF DIPOLE OVER INFINITE PEC GROUND 

The Matlab program to compute the impedance of a dipole in the presence of an 
infinite PEC ground is written based on Equations (II.6) to (II.9) in Chapter II. The 
program is as follows: 


clear 

elf 

% geometry flag: 

% 0 = horizontal, height = ht 
% 1 = vertical, height = ht 
flag=0; 

if flag==0, disp(’horizontal dipole over gp'),end 

if flag==l, disp('vertical dipole over gp'),end 

% total length L in wavelengths 

L=0.45; 

k=2*pi; 

it=0; 

efp=377/4/pi; eep=efp/2; 

% horizontal dipole is side by side case with d=2*h 
if flag==0 
for ht=0.05:.01:l 
it=it+l; 

H(it)=ht; 

d=2*ht; 

uO=k*d; 

ul=k*(sqrt(d A 2+L A 2)+L); 

u2=k*(sqrt(d A 2+L A 2)-L); 

R(it)=efp*(2*Ci(uO)-Ci(ul)-Ci(u2)); 

X( it)=-efp * (2 * S i(uO) - S i(u 1)- S i(u2)); 
end 
end 

% vertical is colinear case with spacing h=(L+s)and d=0) 
if flag==l 

for ht=L/2+.01:.01:1.225 
it=it+l; 

H(it)=ht-L/2; %height of dipole base above gp 
s=ht-L/2; 
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h=L+2*s; 
vO=k*h; 
vl=2*k*(h+L); 
v2=2*k*(h-L); 
v3=(h A 2 -L A 2)/h A 2; 

R(it)=-eep*cos(vO)*(-2*Ci(2*vO)+Ci(v2)+Ci(vl)-log(v3))... 

+eep*sin(vO)*(2*Si(2*vO)-Si(v2)-Si(vl)); 

X(it)=-eep*cos(vO)*(2*Si(2*vO)-Si(v2)-Si(vl))... 

+eep*sin(vO)*(2*Ci(2*vO)-Ci(v2)-Ci(vl)-log(v3)); 

end 

end 

figure(l) 
hold on 

plot(H,R,'k-',H,X,'k-.'),grid 
xlabel('Height above ground plane (wavelengths)') 
ylabel('Resistance or Reactance (ohms)') 
legend('Resistance','Reactance') 

title(['Black: Theory Blue: Patch Red: NEC-Win Pro’]) 


B. ARRAY DIRECTIVITY 

Equations (11.38) and (11.39) define the array directivity derived by classical 
approach. The Matlab program to plot the array directivity as a function of scan angle is 
as follows: 


% This version: array of z-dipoles along x 
clear 

rad=pi/180; 

ntl=3; % number of dipole elements 

h=0.25; % height above ground plane at y=0 

wdth=0.005; % strip width (not used) 

dx=0.5; % separation between dipole elements 

x=[-l:l]*dx; 

y=zeros(size(x)); 

freq=300e6; 

wave=3e8/freq; 

bk=2*pi/wave; 

% dipoles are along the z axis 

% y axis normal to aperture; infinite ground plane at y=0 
%thetas=input('enter scan theta:'); 
thetas=90 
for count=l:91 
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phis=count-1 

us=sin(thetas*rad)*cos(phis*rad); 
vs=sin(thetas * rad) * sin(phis * rad); 

% read integration constants (gaussian quadrature) 
load gausq20.m 
xt=gausq20(:,l); 
at=gausq20(:,2); 
nt=length(xt); 

% set the number of integration intervals (nt points per interval) 
ndivt=2; 
ndivp=8; 

% integration interval in theta (degrees) 

Sl=0*rad; 

S2=180*rad; 

% integration interval in phi (degrees) 

Ql=0*rad; 

Q2=180*rad; 

% generate integration points in theta and phi 
ds=(S2-Sl)/ndivt; 
for i=l:ndivt+l 
SS(i)=(i-l)*ds; 

disp(['i,SS(i)= ’,num2str(i),', ',num2str(SS(i))]) 
end 

dq=(Q2-Ql)/ndivp; 
for i=l:ndivp+l 
QQ(i)=(i-l)*dq; 

disp(['i,QQ(i)= ’,num2str(i),', ’,num2str(QQ(i))]) 
end 

% subintervals in phi 
nphi=0; 
for ii=l:ndivp 
Pl=dq/2; 

P2=(QQ(ii+1 )+QQ(ii))/2; 
for n=l:nt 
nphi=nphi+1; 
wphi(nphi)=at(n); 
phi(nphi)=P 1 *xt(n)+P2; 
end 
end 

% subintervals in theta 
ntheta=0; 
for ii=l:ndivt 
Tl=ds/2; 

T2=(SS(ii+l)+SS(ii))/2; 
for i=l:nt 
ntheta=ntheta+1; 
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wtheta(ntheta)=at(i); 
theta(ntheta)=Tl *xt(i)+T2; 
end 
end 

ninteg=ntheta * nphi; 

disp(['number of int points= ',num2str(ninteg)]) 

% compute field at the integration points 
emax=0; 
sumx=0; 
phase=0; 
for iphi=l:nphi 

% disp( ['phi- ,num2str(phi(iphi)/rad) ]) 
for itheta=l :ntheta 
thr=theta(itheta); 
st=sin(thr); ct=cos(thr); 
phr=phi(iphi); 
cp=cos(phr); sp=sin(phr); 
sumAF=0; 

% direction cosines in the global system 
u=st*cp; v=st*sp; w=ct; 

% element factor sin(theta) 

EF=abs(st); 
for n=l:ntl 

argx=bk*x(n) * (u-us); 
argy=bk*y(n)*(v-vs); 
sumAF=sumAF+exp(j *(argx+argy)); 
end 

% have calculated the E-theta component 

Etheta=sumAF*EF*2*j*sin(bk*h*sin(phr)); 

Ephi=0; 

% all three of these give the same result (as expected) 
emagsq=abs(Etheta) A 2; 
xint=emagsq*sin(theta(itheta)); 

if emagsq>emax, emax=emagsq; thmax=thr/rad; phmax=phr/rad; end 
sumx=sumx+wphi( iphi) * wtheta(itheta) * xint; 
end % end of theta loop 
end % end of phi loop 
prad=T 1 *P 1 *sumx; 
gain=4*pi*emax/prad; 
gdb(count)=10*logl0(gain); 
end 
gdb; 

scan=0:l:90; 

gdb=fliplr(gdb) 

plot(scan,gdb) 
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xlabel('Scan angle (degrees)’) 
ylabel(’Directivity (dB)’) 

grid on 
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